Preprint typeset in JHEP style. - PAPER VERSION 



Chiral field theories, 
matrix models 



HU-EP-03/29 
IFT-UAM/CSIC-03-16 
UCB-PTH-03/13; LBNL-53386 

Konishi anomalies and 



m 

^ ! K. Landsteiner 

^SI I Instituto de Fisica Teorica C-XVI 

Universidad Autdnoma de Madrid 
28O49 Madrid, Spain 
Karl.Landsteiner@uam. es 

C. I. Lazaroiu 

^ ■ Humboldt Universitdt zu Berlin 

Newtonstrasse 15, 124-89 Berlin- Adlershof 
Germany 

calin@physik.hu-berlin.de 



(N 
00 



O 
m 
o 



Oh! 



Radu Tatar 

Theoretical Physics Group 

Lawrence Berkeley National Laboratory 

Berkeley, CA 94720, USA 



^ . rtatar@socrates.berkeley.edu 

> : 
•i— I , 

^ ! Abstract: We study a chiral J\f = I, U{N) field theory in the context of the 

5^ i Dijkgraaf-Vafa correspondence. Our model contains one adjoint, one conjugate sym- 

metric and one antisymmetric chiral multiplet, as well as eight fundamentals. We 
compute the generalized Konishi anomalies and compare the chiral ring relations 
they induce with the loop equations of the (intrinsically holomorphic) matrix model 
defined by the tree-level superpotential of the field theory. Surprisingly, we find that 
the matrix model is well-defined only if the number of flavors equals two! Despite 
this mismatch, we show that the expansion of the loop equations agrees with 
the generalized Konishi constraints. This indicates that the matrix model — gauge 
theory correspondence should generally be modified when applied to theories with 
net chirality. We also show that this chiral theory produces the same gaugino su- 
perpotential as a nonchiral SO{N) model with a single symmetric multiplet and a 
polynomial superpotential. 
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1. Introduction 



A surprising feature oiAf — 1 strong coupling dynamics was uncovered in the seminal 
work of Dijkgraaf and Vafa [1, 2, 3], who found a relation between the gaugino 
superpotential of a confining Af — 1 theory and certain holomorphic [4] matrix 
models. The recipe they proposed takes the tree- level superpotential of such a theory 
to be the action of the dual matrix model. This conjecture was proved for a few 
nontrivial examples, via two distinct methods. One approach [5] uses covariant 
superfield techniques in perturbation theory to integrate out massive matter fields 
in a gaugino background. A different method was proposed in [6, 7], where it was 
shown that the loop equations of the matrix model coincide formally with chiral ring 
relations induced by certain generalizations of the Konishi anomaly. 

Up to now the matrix model - field theory correspondence has been applied 
almost exclusively to the non-chiral case ^ [l]-[27]. The motivation of the present 
paper is to test the conjecture for the case of chiral models. Ideally, one would like 
to know if matrix models can be used to calculate effective superpotentials of SUSY- 
GUTs or other supersymmetric extensions of the standard model. Of course we are 
far from answering this question. Instead, we will study a model with gauge group 
U{N) and chiral matter content chosen to allow for a straightforward large N limit. 
The matter consists of a field $ in the adjoint representation, two fields A, S in the 
antisymmetric and conjugate symmetric two-tensor representations, and eight fields 
Qi . . .Qs in the fundamental representation to cancel the chiral anomaly. The tree 
level superpotential has the form: 

8 

Wtree^tr[W(^)+S^A]+J2Q'fSQf , (l-l) 

/=1 

where W is a complex polynomial. 

This model has the advantage that the number of matter fields is independent 

of the rank of the gauge group, thus allowing for a large N limit with fixed matter 
content. Further motivation to study this model is provided by its interesting type 
llA/M-theory realization [30]-[34] ^. By taking the strong coupling limit (which 
amounts to lifting the brane configuration to M-theory), it was argued in these 
references that the model is described by a smooth curve. This can be interpreted 
as the existence of a mass gap and confinement, a conclusion which is of course also 
suggested by our model's one loop beta function. The geometric engineering of such 
models is discussed in [36] by using methods of [25], [3 7] -[43]. 

^See however [28, 29] for calculations of the effective superpotential using Konishi anomaly 
relations in certain chiral field theories. 

^We note that there also exist JIB brane configurations describing chiral models, for example 
the so-called "brane boxes" of [35], whose stability is unfortunately unclear. Geometric duals for 
these configurations are not currently known. 
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In the present paper, we study the gaugino superpotential obtained after confine- 
ment. We shall show that the effective superpotential agrees with that of a different 
M = 1 theory, namely a non-chiral SO{N) model with a single chiral superfield X 
in the symmetric representation and a tree-level superpotential given by iiW{X) 
(the gaugino superpotential for such models was recently investigated in [18]). This 
relation can be understood most easily by turning on a D-term deformation of the 
original model, under which the theory flows at low energies to the SO{N) model 
with symmetric matter. Because the effective superpotential is protected by holo- 
morphy, its form must be independent of the choice of Fayet-Iliopoulos term, which 
explains why one obtains agreement between the two theories. One can make this 
argument more precise by computing the first (a.k.a. Veneziano-Yankielowicz) ap- 
proximation to the gaugino superpotentials upon using scale matching techniques, 
and we shall do so below, finding agreement. To give a complete proof of low-energy 
equivalence, we use the more powerful technique of chiral ring relations [7], which 
allows us to characterize the exact effective superpotential in terms of solutions to 
certain algebro-differential equations induced by generalized Konishi anomaly con- 
straints. Then the connection between the chiral and SO{N) models follows upon 
matching the relevant relations in the chiral rings. 

In principle, the generalized Konishi constraints derived below suffice to com- 
pletely determine the exact gaugino superpotential, which can be extracted with 
arbitrary precision by solving the relevant equations. However, it is interesting to 
follow the beautiful insight of Dijkgraaf and Vafa in order to construct a holomorphic 
matrix model whose free energy specifies the superpotential. This can be achieved by 
building a matrix integral whose loop equations reproduce the Konishi constraints. 
Applying these ideas, one finds some novel phenomena, which are related to the 
chiral character of our matter representation. 

In fact, the holomorphic matrix model with action given by the tree level super- 
potential of our chiral field theory turns out to be ill defined. The problem is that, 
although the matrix model action is invariant under the complexified GL{N, C) gauge 
group, the measure fails to be invariant unless the number of matrix model flavors 
Qf equals two. This phenomenon, which is due to the presence of a chiral matter 
content, forces us to work with a matrix model which contains only two flavors, even 
though the associated held theory contains eight ! Then the matrix partition func- 
tion is well-defined, and we show that the loop equations agree with the generalized 
Konishi constraints despite the mismatch in the number of fiavors. Our matrix model 
is intrinsically holomorphic, in the sense that it docs not admit a real or Hermitian 
version. This is due to the fact that our matter representation is chiral. 

Having extracted the relevant matrix model, we compare it with the model which 
governs the gaugino superpotential of the SO{N) theory with symmetric matter, a 
comparison which sheds different light on the relation between the two theories. 

The paper is organized as follows. In Section 2, we analyze the classical moduh 
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space of our theory by solving the F- and D-flatness constraints. For a diagonal vev 
of $, we find that the gauge group is broken to a product YYi=i U{Ni), where each 
factor contains the same massless matter as the original theory. Computation of 
the leading contribution to the effective superpotential requires threshold matching, 
which cannot be performed directly since chirality forbids the addition of a mass 
term. This problem was also encountered in [28], where it was solved by deforming 
the superpotential in such a way as to Higgs the gauge group. We will use a similar 
technique. Instead of deforming the tree-level superpotential, use independence of 
the effective superpotential of D-term deformations, which allows us to add a Fayet- 
Iliopoulos term. In the presence of a Fayet-Iliopoulos parameter (which we take 
to be positive), we find that the symmetric field acquires a vev, which breaks the 
gauge group down to SO{N) with massless matter X transforming in the symmetric 
representation and a superpotential trW^(X). Then the first approximation to the 
gaugino superpotential can be computed by standard threshold matching. The re- 
sulting nonchiral SO{N) model was recently studied in [18] via the Konishi anomaly 
approach of [7]. 

In Section 3, we derive the generalized Konishi constraints for our chiral theory 
and show how they relate to those extracted in [18] for the SO{N) model. Since 
our chiral model has a different matter content, we find a different set of resolvent- 
like objects which enter the relevant 'loop equations'. However, we show that all 
such quantities are uniquely determined by the solution of a pair of equations which 
coincide with those derived in [18] for the SO{N) theory with symmetric matter. 

Section 4 discusses the matrix model dual to our chiral theory. We show that 
the measure is not invariant under the central C* of the GL{N,C) gauge group (and 
thus the matrix model partition function vanishes or is infinite) unless the number 
Np of matrix model fiavors equals two. We then extract the loop equations of this 
model by using both the standard method of the eigenvalue representation and the 
approach of [22]. Finally, we discuss the relation with the matrix integral relevant 
for the SO{N) theory with symmetric matter. 

The identifications mapping our loop equations into the Konishi constraints are 
given in section 5. Using this map, we extract an explicit formula expressing the 
gaugino superpotential in terms of the matrix model free energy. This completes the 
proof of the Dijkgraaf-Vafa conjecture for our case. 

Section 6 presents our conclusions. In appendix A we prove gauge-invariancc of 
the matrix model measure. Appendix B recalls the classical vacua of the SO{N) 
theory with symmetric matter, while appendix C contains some details relevant for 
the discussion of Section 3. 

2. A first view of field theory properties 

In this section we take a first look at our field theory model. After describing it 
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precisely, we discuss the part of the classical moduli space which will be relevant for 
our purpose, and give our derivation of the Veneziano-Yankielowicz superpotential, 
which is the first approximation to the exact glueball superpotential predicted by 
the Dijkgraaf-Vafa correspondence. 

2.1 Description of the model 

We start with a U {N) gauge group, together with chiral matter $, S, A in the adjoint, 
antisymmetric and conjugate symmetric representations, as well as Np quarks Q/ in 
the fundamental representation. We consider the tree-level superpotential: 

Np 

Wtree^tT[W{<^)+S<^A]+Y.QPQf ' (^-l) 

/=1 

where: 

d+l 

W{z) = J2 (2-2) 
i=i 

is a complex polynomial of degree d + l. We have S'^ — SjA^ — —A (while $ is 
unconstrained) and the gauge transformations are: 

^^C/dC/^ , S^USU^ , A^UAU^ , Qf^UQf , (2.3) 

where U is valued in U{N). The U{N) gauge symmetry is obviously preserved by 
Wtree- Notc that the fields S, A are complex. 

We do not have quarks in the anti-fundamental representation and therefore 
this system is quite different from models studied in [8, 14, 15, 16]. In fact, the 
matter representation is chiral, in contrast to most situations previously studied in 
the context of the Dijkgraaf-Vafa conjecture. In particular, the model will have a 
chiral anomaly unless we take Np — 8. In the following, we shall focus on the non- 
anomalous case though we allow Np to take an arbitrary value in most formulas (this 
permits us to recover the anomaly cancellation constraint Np = 8 as a consistency 
condition required by gaugino condensation). 

This model can be obtained through an oricntifoldcd Hanany-Witten construc- 
tion [30, 31, 32]. It can also be realized through geometric engineering, as we discuss 
in a companion paper [36]. 

2.2 The classical moduli space 

Let us study the classical moduli space of our theories. Part of the discussion below 
is reminiscent of that given in [37, 42] and [21] for quiver gauge theories, though of 
course we have a rather different matter content and hence the details are not the 
same. 
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The F-flatness constraints are: 

^^S^S^ , (2.4) 

Np 

^A-A^^ + 2j2QfQ^f^0 , (2.5) 

/=i 

W'{^) + AS^O , (2.6) 

SQf^O , (2.7) 

while the D-flatness condition is: 

1 1 

^[^\^] + S^S-AA^--J2QfQ\-0 , (2.8) 

/=i 

where the left hand side is the moment map for our representation of U{N). 
To understand the solutions, notice that (2.6) implies: 

W'i^f = {ASf = AS AS . (2.9) 

Using the transpose SA — W'{^'^) of (2.6) in the right hand side gives: 

Vr'($)2 = AW'{<I>^)S = ASW'{<I>) (2.10) 

where in the last equality we used equation (2.4). Applying (2.6) once again in the 
right hand side of (2.10), we find: 

W'i^y^O . (2.11) 

Let us assume that [$^, $] = 0, i.e. $ is a normal matrix. Then $ is diagonal- 
izable via a unitary gauge transformation, and equation (2.11) shows that $ can be 
brought to the form: 

$ = diag(Ail^,...Adl^J (2.12) 

where Ai . . . are the distinct roots of W'{z) and Ni . . . and non-negative integers 

such that Ni-\ \- N^ — N. li Nj — for some root A^, we use the convention that 

the corresponding block XjInj does not appear in (2.12). 

With this form of equation (2.4) shows that 5* must be block-diagonal: 

S^disigiS,...Sd) (2.13) 

where Sj are symmetric Nj x Nj matrices. When bringing $ to the form (2.12), we 
are left with a residual 11^=1 ^i^j) g^^^gc symmetry corresponding to the transfor- 
mations U — diag(C/i . . .Ud) with Uj G U{Nj). Using this symmetry, we can bring 
Sj to the form^: 

Sj = diag(0^(o), <7]^)l^a) . . . (2.14) 



^This is the Takagi factorization (see, for example, [44] page 204) for the complex symmetric 
matrix Sj. 
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where: 







< (7, 



(rrij) 



(2.15) 



The multipUcities ivj'^-' satisfy Ylh^o 



Ni. In the case m,- = 0, we have N) 



(0) 



Nj and equation (2.14) reduces to Sj = 0^^. After bringing S to the form (2.14), we 
left with 
Writing: 



are left with the gauge symmetry if U{Nf) nr^i SO{N^''^) 



Qf 



(1) 



Qf 



(2.16) 



where Q^p is an A'j-vector, equation (2.7) requires that each lies in the kernel 



of Sj. Thus we must have: 



Qf 



rJJ) 







(2.17) 



when decomposing into sub- vectors according to Nj — N^^^ + Nj^'' + ■ ■ ■ + 
Here qf is a column iVj°^-vector. Using the remaining gauge symmetry, we bring 



the X Np matrix q^^ :-- 



Uf 



S3) 



to the form: 



.0) .- 




* * 



a,. * ... * 
... 



(2.18) 



... 

where < ai < • • • < a^^ , the symbol * stands for generally distinct complex entries 
and the zero rows at the bottom may be absent. The rank Sj of this matrix equals the 



dimension of the vector space spanned by q\ 



and mm{N^^\ Np). If we let pj := - Sj, then the vevs (2.18) break the C/(A/"f 
components of the gauge group down to U{pj). 

Finally, equation (2.5) shows that in such a vacuum the x Nj blocks of A are 
given by: 

2 



r(0) 



q^^ , and can take values between 

{0)^ 



A. 



3 « /=1 



(2.19) 



Since Qf^ lie in the kernel of Si., this automatically satisfies condition (2.6), which 
in our vacuum takes the form AS — <(=^ SA — 0. Equation (2.19) shows that 
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all non- vanishing entries of A lie in the Sj x sj sub-blocks of the Nj^'^^ x Nj'^^ blocks, 
where they are given by equations (2.19) and (2.18). 

Since we assume [$^, $] = 0, the D-fiatness condition (2.8) reduces to: 



S^S-AA^ --Y,QfQ\-0 ■ (2.20) 
/=i 

Using the form of $, S and A discussed above (which is required by F-flatness), it is 
not hard to see that (2.8) implies (note that we do not sum over i): 

SjSi^O foralH . (2.21) 

This follows by decomposing (2.20) into iVj x Nj blocks and restricting to the case 
i = j while using the fact that {AA'^)ij — —{AA)ij = —Ai^Akj and (Q/Qy)?^ = 
both vanish for i = j due to equations (2.16) and (2.19) except in the A^^^^^ x A^^^'^^ 
block. This implies that the entries of the block diagonal matrix SiSj vanish except 
in this block. The entries also vanish in the N^^^ x N^^^ block because of the form 
(2.14). 

Since SjSi is positive scmidefinite, equations (2.21) imply Si = for all i and 
thus S — 0. In this case, the D-flatness condition (2.20) becomes: 

AA^ + -J2QfQ\ = , (2.22) 
/=i 

which implies A = and Qj = for all / by semi-positivity of the left hand side. 

It follows that the only classical vacua for which [$^^, $] = are given by $ of 
the form (2.12) and S = A = as well as Q/ = for all /. In such a vacuum, the 
gauge group is broken down to the product 11^=1 ^i^j)- 

2.3 The Veneziano-Yankielowicz superpotential 

We next discuss the leading approximation to the gaugino superpotential. As we 
shall see below, the effective superpotential coincides with that of an SO{N) field 
theory with a single chiral superfield X transforming in the symmetric two-tensor 
representation, and a tree-level superpotential trW^X). 

To compute the gaugino superpotential, we need the scalc(s) of the low en- 
ergy theory, which are usually obtained via threshold matching. Standard threshold 
matching is difficult to apply to chiral theories, since a chiral tree-level action cannot 
contain mass terms for the chiral fields, and thus one cannot directly integrate out 
such fields at one-loop. This problem was encountered for a chiral model considered 
in [28], where it was overcome by applying threshold matching to a certain Higgs 
branch. Some aspects of the same issue were recently discussed in [19]. 
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It turns out that threshold matching can be carried out in our case provided 
that one first deforms the theory through the addition of a Fayet-Ihopoulos term ^ 
(this is useful for our purpose since holomorphy dictates that the low energy scale is 
insensitive to D-term deformations) . Such a deformation was previously considered in 
[30, 31, 32] and has the following effect. Concentrating on vacua with [($), ($)''"] = 0, 
it is not hard to see that the inhomogeneous form of (2.8) (obtained by introducing 
^Iat in the right hand side) together with the F- flatness constraints again imply 
(A) — and (Qf) = for all /. On the other hand, the symmetric field S gets 
an expectation value equal to the square root of ^, which breaks U{N) to SO{N). 
The vev of S gives equal masses to the quarks Qf due to the terms Q^SQf. The 
term tr S^A gives equal masses to the fields A and Y where F = |(<l> — is the 
antisymmetric part of $. All of these masses depend on the Fl-term. The symmetric 
part X := |($ + $^) of the adjoint field remains massless. The fluctuations of 
S around its vev are 'eaten up' while giving masses to the appropriate W-bosons 
through the Higgs mechanism. 

Eventually, one is left with an SO{N) theory with a symmetric tensor X and 
a supcrpotential tTW{X). By holomorphy, the scale of this SO{N) theory must be 
independent of the original Fl-parameter. It is easy to check this explicitly. Let 
the Fl-paramctcr be = for some real v. Then the inhomogeneous form of (2.8) 
shows that the vev of S equals ±z/, and we can take the plus sign without loss of 
generality. The massive fields are eight quark fiavors, two antisymmetric tensors (A 
and the antisymmetric part of $) and the W-bosons in the symmetric representation 
of the low energy gauge group SO{N). This gives the scale matching relation: 

A^-4 = A^-V4i/^-2i/-(^+2) = A^-4 ^ Ao = A . (2.23) 

Here A is the scale of the chiral high energy theory (whose one-loop beta function 
coefficient is — 4) and Aq is the scale of the low energy SO{N) theory. Notice 
that the exponent of Aq is unusual in the sense that we would expect 2A^ — 8 for 
the SO{N) theory with symmetric tensor. This can be traced back to the fact that 
the generators in the SO{N) theory are unusually normalized. It can be seen by 
noticing that the index for the fundamental representation of U {N) is | whereas the 
index for the fundamental representation of SO{N) with conventional normalization 
is 1 and that a fundamental of the high energy U (N) theory descends directly to a 
fundamental of the low energy SO{N) theory. This normalization has already been 
taken into account in (2.23). 

We still have the deformation tYW{X), which leads to diagonal vevs of X of 
the form (2.12). The relevant vacua of the SO{N) theory with a symmetric field are 
discussed in Appendix B. One finds that the vev of X further breaks the Lie algebra 
of the low energy gauge group according to [45] : 

so{N) ^ ®'^^^so{Ni) . (2.24) 
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It is now easy to extract the scales of the different so{Ni) factors (here we use 
conventional normalization, since we compare only SO theories): 



mv' I I m^,.. 



n 



A2(^-4) (2.25) 



where ruxi = W"{Xi) and the masses of the SO{N)/ SO{Ni) W-bosons are mvy.^ = 
Xi — Xj. Because of (2.23), wc can use the high energy scale A in (2.25). 

The Veneziano-Yankielowicz contribution to effective superpotential has the form: 

d <^ /.3(7Vi-2)\ 

Weff = E y [log ) + iV. - 2] , (2.26) 



=1 



where we have inserted the factor 1/2 coming from the relative normalization of the 
generators of U{N) and SO{N). 

Observation: One can also consider turning on a negative Fl-term ^ = —v^. Tak- 
ing to be even for simplicity and assuming [($)^ ($)] = and (Q/) = for all /, 
the D-flatness condition shows that the antisymmetric field A acquires a vev. Up to a 
gauge transformation, we can take (A) = vJ where J is the antisymmetric invariant 
tensor of ,5p(A^/2): 

This breaks the C/(A") gauge group down to Sp{N/2), with the fluctuations of A 
'eaten' by the W-bosons of the coset U{N) /Sp{N/2). Let us decompose $ = 
where ($±J)^ = ±($±J). Then the superpotential becomes: 

tr[W^($+ + <l>_) + z/5$+J + Q^®gJS] . (2.28) 

Integrating out $_|_ by the equation of motion of S (which gives <l>_|_ = ~j;Qf ® Q^J) 
results in the superpotential: 

Wio^ = W{^.-^Qf0Qp) . (2.29) 

Hence at low energies we have an Sp{N/2) gauge theory with an antisymmetric tensor 
and eight fundamentals interacting through the superpotential (2.29). Because of 
the complicated structure of Wiow, we did not find this branch to be useful for our 
purpose. 

3. Low energy analysis via generalized Konishi anomalies 

In this section, we extract the relevant chiral ring relations of our model and compare 
with those of the SO{N) theory with symmetric matter. We shall use the method 



10 



of generalized Konishi anomalies originally developed in [6, 7]. The structure of the 
tree level superpotential implies: 



dW ft 

J^ = (tr($^)) . (3.1) 

Our strategy is to extract a set of Konishi anomaly relations which allow one to 
solve for the generating function T{z) — (tr {jz^)) of the chiral correlators (tr 
appearing in the right hand side. The integration of (3.1) allows one to compute the 
effective superpotential up to a piece which is independent of the coupling constants 
tj- 

3.1 Konishi constraints for the chiral model 

As discussed in [7], the loop equations of the (adjoint) one-matrix model are formally 
equivalent to certain chiral ring relations induced by a generalized form of the Konishi 
anomaly of the U{N) field theory with one adjoint multiplet. The argument extends 
to other matter representations, and is based on special properties of chiral operators 
(^operators corresponding to the lowest component of chiral superfields) in jV = 1 
supersymmetric field theories in four dimensions. It is well-known that correlators of 
such operators do not depend on the space-time coordinates and that they factorize: 

(OiO2) = (0i)(02) . (3.2) 

On the set of chiral operators one considers the equivalence relation: 

Oi = 02 + Cag"(...) , (3.3) 

where c^Q" is an arbitrary linear combination of the anti-chiral supercharges. This 
equivalence relation is compatible with the operator product structure, and modding 
out the space of chiral operators by (3.3) leads to the chiral ring. Equivalence of two 
chiral operators under (3.3) implies equality of their vevs: 

m = {0,) . (3.4) 

An important relation in the chiral ring is: 

{>V(r),<^}^0 , (3.5) 

which holds for any representation r of the gauge group. This is a special case of the 
general relation [28] : 

yV;W.c)W = o , (3.6) 

which holds for an arbitrary chiral operator O^'^^ transforming in the representation 
r of the gauge group. The dot in this equation indicates the action of Wq in the 
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representation r. For our theory, one finds: 



= m,^=0 (3.7) 

W^.A = WaA + AWj = (3.8) 

Wa.S = -SWa - W^S = (3.9) 

Wa-Qf ^yVaQf = , (3.10) 



where Wq in the right hand sides are taken in the adjoint representation and jux- 
taposition stands for matrix multiphcation. Note that Wq, acts on the product AS 
through the commutator (since this product transforms in the adjoint representa- 
tion) . 

The generalized Konishi anomaly is the anomalous Ward identity for a local 
holomorphic field transformation: 



The supercurrent generator has the form J = (0^^))'''e^*'^'50*^'') where V^^^ is the 
vector superfield in the representation r and O*^'') is the chiral superfield associated 
with The chiral ring relation induced by the generalized Konishi anomaly for 

this current is: 

where the capital indices enumerate a basis of the representation r. We will investi- 
gate the generalized Konishi relations corresponding to the field transformations: 

5$ = J^lJZ^ (3.13) 

5$ = (3.14) 

M ^ (3.15) 

= -^"^—^ (3-16) 
z - $ z - $^ 

= —^^^^ (3-17) 

Np 

^Qf = E JZ^Q^ ■ (3-18) 

9=1 



In the last equation, A is an arbitrary matrix in flavor space. 
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Writing = WWa, we define: 

Then it is shown in Appendix C that transformations (3.13-3.18) generate the chiral 
ring relations: 

: '1 - T^T^tr S-^A ^ R{zr (3.22) 



tr (^^^-^] +tr f ) = 2R(^)T(z) + w"(^)w«(^) (3.23) 



z - ^ J \ z - ^ 

1 \ 1 



,tr ( S^^A ] = -Rizf (3.24) 
327r2 \ z-<^ J 2 ^ ^ ^ ^ 

tr (^Sj^A^ = R{z)T{z) + 2R'{z) - ^w"(^)w„(^) (3.25) 
i^J^^) + ^ ^^^^^^^^ ^ R'(^)T(^) - 2R'(^) - ^w'^(^)w.(^) (3.26) 

Taking the trace of the last equation gives: 

2Q^S-^Qf = R{z)NF . (3.28) 

We next take the vacuum expectation values of these chiral ring relations. Let 
us define: 

R(z) := {R(z)) (3.29) 
T{z) := {T{z)) (3.30) 

M{z) := (tr (^S-^A^^^ (3.31) 
Mq{z) = J2 {^^fJ^^J^^f) (3.32) 

^(^)^=-3^f(^i^^)) ^'-''^ 
L(^) := J] /gj5^g A . (3.34) 
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Introducing the degree d — 1 polynomials: 

= i ^ ' (3-36) 



we write: 



327r2 



1 ^ - /(^) (3.37) 

(yrf ) = W^'(^)n^) - c{z) . (3.38) 

Noticing that the vevs of spinor fields vanish due to Lorentz invariance, we find 
the following Ward identities for the generating functions (3.29-3.34): 

R{z)^ - K{z) - W\z)R{z) + f{z) = (3.39) 
2R{z)T{z) - W'{z)T{z) - M{z) + c{z) = (3.40) 

K{z) - \R{zf = (3.41) 

M{z) - R{z)T{z) - 2R'{z) = (3.42) 
Mq{z) + M{z) - R{z)T{z) + 2R'{z) = (3.43) 
2L{z) - NfR{z) = . (3.44) 

It is easy to eliminate K{z), M{z), Mq{z) and L[z) from these equations to find: 

R{zf - 2W'{z)R{z) + 2f{z) = (3.45) 
T{z)R{z) - W'{z)T{z) - 2R!{z) + c{z) = . (3.46) 

Given a solution {R{z),T[z)) of these constraints, the quantities K,M,Mq and L 
are given by: 

K{z) = ^Rizf 

M{z) = 2R!{z) + R{z)T{z) (3.47) 

Mq{z) = -m'{z) 

Np 

L{z) = ^R{z) . 

Hence all solutions are parameterized by the 2d complex coefficients of the polyno- 
mials f{z) and c{z). 

The generalized Konishi relations involving the fiavors Qf have an interesting 
implication. Expanding the last two equations in (3.47) to leading order in 1/z 
gives: 

J^iQpQf) = 4<5 (3.48) 
/ 
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and 

Y.^QpQi) = , (3.49) 

/ 

where S = — 3^(trW^) is the gaugino condensate. If S is non- vanishing, then 
compatibihty of these two equations requires that we set Np — 8, which is also 
required for cancehng the chiral anomaly. Any other value is incompatible with the 
existence of a gaugino condensate. 

Observation: The particular form of (3.27) is related to the 0{Np) flavor symme- 
try of the field theory with tree-level supcrpotential (2.1): 

Qf^Q'f-.^TfgQg , (3.50) 

where r is a general complex orthogonal matrix. Since this symmetry is Tinbroken 
after confinement, the quantity Lfg = (Q^S^^Qg) must be 0{Nf) invariant and 
thus proportional to Sfg. Equation (3.27) shows that the proportionality factor is 
given by ^R{z) . 

3.2 Compcirison with the SO{N) model with symmetric matter 

The Konishi relations (3.45) and (3.46) coincide with those of the SO{N) field theory 
with a single complex chiral superfield X transforming in the symmetric representa- 
tion, and with a tree-level superpotential given by trW{X). The Konishi relations 
for this theory were derived in [18] (see equation (18) of that reference^) and one 
immediately checks that they agree with our equations (3.45) and (3.46). More pre- 
cisely, the SO{N) theory with symmetric matter admits the following resolvent-like 
objects: 

Rs{z):^-{-^tr (^]) (3.51) 



Tsiz) := (tr (j^)) , (3.52) 

which were shown in [18] to obey the two equations (3.45) and (3.46). This leads to 
the identification: 

R{z) « Rsiz) 

T{z) ^ T,{z) . (3.53) 

which maps a solution of our Konishi constraints to a solution of the Konishi relations 
of the SO{N) theory (this works because equations (3.47) completely determine the 
other resolvent-like objects of the chiral theory given a solution {R, T) of (3.45) and 
(3-46) ). 

^The paper [18] considers a slight extension of our SO{N) theory, by adding some fundamental 
and anti-fundamental matter beyond the symmetric field. The case of interest for our purpose is 
recovered from the equations of [18] by setting their quark mass matrix to zero. 
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4. The matrix model 



The general conjecture of [1] suggests that the effective sup ci potential of our field 
theory should be described by the holomorphic matrix model^: 

Z=-^ [ d^dAdSdQe-^^"'"'^^'^'^'^^ , (4.1) 
1^1 Jr 

where IGI is a normalization factor and: 



A, S, Q) = tr \W{^) + S^A^ +Y, QpQf = tr 

/=i 



Np 

/=i 

Here $ is an arbitrary complex N x N matrix, the complex N x N matrices S 
and A are symmetric and antisymmetric respectively and Q is a general complex 
N X Np matrix whose columns we denote by Q/- The integration measure d/x = 

d^dAdSdQ = d^dAdSYl^Zi dQf is given by: 

N Np N 

d^= /\d^ij , dS = f\dSij , dA^/\dAij , dQ^/\/\dQ} . (4.3) 

i,j=l i<j i<j f=l 1=1 

where /\ denotes the wedge product and we use the lexicographic order of indices to 
give unambiguous meaning to the various products of one-forms. For example, the 
notation Ai<j ^.Sij means: 

dSiiAdSuA- ■ ■ AdS^j^ AdS22 /\dS23A- ■ ■ AdS2ffA- ■ ■ AdSj^_^j^_-^^AdSji;f_^j^ AdS]<^ji;f . 

(4.4) 

Of course, the ordering convention can be chosen arbitrarily since changing it pro- 
duces an irrelevant sign pref actor in the matrix integral. 

The measure d^ is a top holomorphic form on the complex space: 

M = {{^,S,A,Q^...Qf,^)\S'' = 3,^" = -A} . (4.5) 

The integral in (4.1) is performed on a boundary-less real submanifold T ol M. whose 
closure is non-compact and which is chosen such that dimigr = dimcA^. The model 
(4.1) admits the 0{Nf) flavor symmetry: 

Qf^Q'r-=rfgQ, , (4.6) 
where r is an Np x Np orthogonal matrix . 



^It is interesting to notice that this holomorphic matrix model does not have a real ('Hermitian') 
counterpart. This is due to the fact that our matter representation is intrinsically complex. 
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Since anomaly cancellation in our field theory requires Np = 8, one is tempted 
to set Np = 8 as well. It turns out that this naive identification cannot hold in our 
case. To understand why, notice that both the matrix model action (4.2) and the 
integration measure are invariant under the following SL{N, C) gauge transforma- 
tions: 

$ ^ Um-^ , S {U-^fSU-^ , A UAU^ , Qf UQf , (4.7) 

where C/ is a complex N x N matrix of unit determinant (invariance of the measure 
is discussed in detail in Appendix A). To preserve this symmetry, one must choose 
r to be stabihzed by the action (4.7) of SL{N, C). 

The matrix model action is in fact invariant under the full GL{N, C) group 
acting as in (4.7). However, the measure dfi is not invariant under the central C* 
subgroup of GL{N,C) unless Np = 2. Taking ^7 = ^1^ in (4.7) with ^ e C*, we 
have: 

A^eA , S^r'S and Qf^^Qf , (4.8) 

which gives: 

^ ^iv(iv^-2)^^ . (4.9) 

Let us assume Np ^ 2 and choose F to be GL(N,C) invariant. Then, since 
the matrix model action (4.2) is invariant under (4.8), while the measure transforms 
nontrivially, invariance of the integral (4.1) under coordinate transformations shows 
that: 

2^^NiNF-2)2 . (4.10) 

Thus Z must either vanish or equal complex infinity^! A similar argument shows that 
the integral J d/iFe^'S'^""^ must vanish or be infinite for any functional F($, A, S, {Qf}) 
which is invariant under (4.8). In particular, the expectation value: 

(F) — ^ / duFe'^^"''^ (4.11) 

of any such functional is ill-defined ! This means that the matrix model predicted 
by a naive application of the conjecture of [1] is not well-defined. 

That subtleties can arise when attempting to apply the conjecture of [1] to chiral 
field theories is not completely unexpected, since most derivations of this conjecture 
up to date have concentrated on real matter representations, which prevent the 
appearance of net chirality. The phenomenon we just discussed shows that one must 
modify the original conjecture of [1] in order to adapt it to the chiral context. 

Thus we arc lead to consider the matrix model with Np = 2. Then both the 
action (4.2) and the integration measure arc invariant under GL{N,C) transforma- 
tions of the form (4.7), where U is now an arbitrary complex invertible matrix. In 

^The second solution is allowed since Z is complex and the point at infinity in the complex plane 
does satisfy oo = ^*(*p-2)oo. 
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Subsections 4.3 and 4.4 below, we shall show explicitly that the model with Np = 2 
is well-defined by relating it to the holomorphic matrix model associated with the 
SO{N) theory with symmetric matter. 



4.1 Loop equations 

In this subsection, we extract the loop equations of the model (4.1). Although the 
correlation functions are not well defined unless Np = 2, we will work formally with 
an arbitrary value of Np. This will allow us to re-discover the constraint Np = 2 as 
a consistency condition between the loop equations, in a manner similar to the way 
in which we recovered the condition Np = 8 in Subsection 3.1. by using the Konishi 
constraints of the field theory. 

In addition to the matrix model resolvent: 



u;{z) = -^tr- 



N z-^ 



(4.12) 



we shall consider the objects: 



k{z) 



N 



-tr 



S- 



2; - $ 



z-(^^ 2; - $ 

l{z) = Q^S-^Q^ . 
z — <v 

We will show that these fulfill the loop equations: 



(4.13) 
(4.14) 
(4.15) 



{u:{zf-W'{z)u{z)^k{z) + f{z)) = Q 

{u'{z)+mQ{z))^Q 
{Npuj{z) - 2l{z)) = , 



(4.16) 

(4.17) 

(4.18) 
(4.19) 



where: 



^ ' N z-^ 



(4.20) 



is a polynomial of degree d — 1. 

Before giving the derivation of these constraints, let us note that one can elimi- 
nate {k{z)) between (4.17) and (4.16) to find an equation for the resolvent: 



uj{zf - -|cu'(^) - 2W'{z)uj{z) + 2f{z] 







(4.21) 
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Given a solution {uo{z))^ relation (4.17), (4.18) and (4.19) determine the averages of 
k{z), mq^z) and l{z) as follows: 



1 
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{k{z)) = -{co{zy + ^u;'{z)) 
{ruQiz)) = -{J{z)) 

{l{z)) = ^{u^iz)) . 

The leading order in the large z expansion of the last two equations gives: 

(gpQ;) = g 

where we used the large z behavior of the resolvent: 

{u;iz)) ^ ^ + Oi-) . 



(4.22) 



(4.23) 
(4.24) 

(4.25) 



Since we of course take g ^ 0, equations (4.23) are consistent only if Np — 2. We 
now proceed to give the proof of (4.16-4.19). 

4.2 Direct derivation of the loop equations 

Consider the identity: 



/ 



Using: 



this leads to: 



d^dAdSdQ 



d 



d 



1 



z - ^ 



cuizf ^ 4tr 

N 



z - ^ 



z - ^ 



-k{z) 







(4.26) 



(4.27) 



(4.28) 



Equation (4.20) allows us to write (4.28) in the form (4.16). 

In order to find an additional equation for k{z), we consider the identity: 



/ 



d^ dA dSdQ 



d 



' -A- ' 



z-^ Z-^'^Jij 



. 



(4.29) 



This implies equation (4.17) upon using the relations: 



$ . 1 
tr I S^^A- 



z - ^ z -^'^ 



1 . 

-tr I S^^A- 



z-^ z - 



tr S 



-A 



(4.30) 
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which follow trivially from the symmetry properties of S and A and invariance of 

the trace under transposition^. 
We next consider the identity: 



d^dAdSdQ 



d 



-s- 



, 



(4.31) 



which gives: 



^u{zr - llu'iz) - k{z) - jmQiz)^ = . (4.32) 



Together with (4.17), this implies the third loop equation (4.18). 

Finally, we can derive a relation involving the flavors. For this, we start with the 
identity: 



/ 



d^dAdSdQ 



d 



dQ 



f 



A 



fg 



Qie 



, 



where A is an arbitrary matrix in flavor space. 
Equation (4.33) gives: 



A/M^)-2^gp-^g,^ = o 



Since A is arbitrary, this implies: 



(4.33) 



(4.34) 



{u;{z)5f,-2Qp ^g,) = 

z - ^ 



(4.35) 



Setting f = g and summing over / gives equation (4.19). 

Since the matrix model admits the 0{Nf) flavor symmetry (4.6), it follows that 
the expectation value of any flavor two-tensor must be 0{Np) invariant and thus: 



1 



1 



z-^ Nf 



(4.36) 



Combining with (4.19), this gives: 



1 



{QjS Q^)= {uj{z))6f, , 

z-^ 2 

which is the matrix model analogue of equation (3.27). 



(4.37) 



^For this, it is useful to notice that tr (S^^A — ^) = 



20 



4.3 The eigenvalue representation 

As mentioned above, the integration in (4.1) should be performed over an appropriate 
multidimensional contour inside the space of complex matrices S, A and Qf. We 
shall specify this contour by first fixing the gauge through diagonalizing ^ $ and 
imposing conditions on the remaining matrices after gauge-fixing. This procedure is 
clearly gauge- invariant in the sense that it defines a multidimensional contour which 
is stabilized by the action of the complex gauge group GL{N, C). Using the gauge 
symmetry to bring $ to the form: 



l' = diag(Ai...A^) 

allows us to write the matrix model action (4.2) as: 



(4.38) 



N 






Nf 


TV 


Nf 










4- + E 




i=l 






/=1 


1=1 


/=1 



To make sense of the remaining integral, we impose the conditions: 
Ai e 7 , Sij e iR , i^^ e K , Q} e K , 



(4.39) 
(4.40) 



where 7 is an open contour in the complex plane whose asymptotic behavior is 
dictated by the highest degree term in as explained in [4] (one can take 7 to 
coincide with the real axis if and only if is a polynomial of even degree). This 
amounts to choosing: 



r = S,A,Qi...Q^^) eM\3U e GL{N,C) such that 
U^U^^ = diag(Ai . . . A^) with Ai . . . Ajy G 7 and 
U-^SU-^ e Mat{N,N,R), UAU^ e Mat{N,N,R), 
UQfe Mat{N, 1, M) for all /} . 



With this choice of integration contour, we obtain: 
where: 

„ N Nf 

Zq= I __dQ\[5{Y,{Q)?) 

j=i /=i 



(4.41) 



(4.42) 



(4.43) 



®The fact that $ is diagonalizablc is part of the definition of our multidimensional contour. 
As in [4], this implements a point-splitting regularization of the matrix integral, which can be 
removed trivially because the resulting Vandermonde determinant (to be described below) vanishes 
for coinciding eigenvalues. 
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and: 



Zred = ld\^...\d\^W {\ - A,)e-f . (4.44) 



To arrive at (4.42), we performed the integrals over Sij^ which appear linearly in the 
action. This gives a product of delta-functions which allows us to reduce the integral 
to (4.42). From this expression, it is clear that the variables A, S and Q/ decouple 
from $. The interesting dynamics of the model is contained in the reduced partition 
function (4.44), which differs from that of a usual (adjoint) one-matrix model only 
because the Vandermonde determinant A = {\ ~ \) is not squared in (4.44). 
As we shall see below, Z^^,^ can in fact be identified with the partition function of a 
holomorphic SO{N, C) - invariant one matrix model with a single symmetric field X 
and action tr W{X). This, of course, is just the holomorphic matrix model associated 
with the SO{N) theory with symmetric matter via the Dijkgraaf-Vafa conjecture. 

Observation: The integral (4.43) is finite and non-vanishing precisely in the case 
of interest Np — 2. In this case, one easily checks that: 

ZQ = 7r^ . (4.45) 

For Np > 2 we find Zq — 0, while for Np — 1 we have Zq — oo. This agrees with 
our previous discussion. 

4.4 Relation to the matrix model of the SO{N) theory 

As mentioned above, it turns out that the reduced model described by (4.44) agrees 
with the matrix model associated with the SO{N) theory with symmetric matter via 
the Dijkgraaf-Vafa correspondence. The Hermitian version of the latter is defined 
through: 

^.:=^/ciXe-f^^^(^) , (4.46) 

where \Gs\ = vol(5'0(iV, M)/5'^) and X is a real symmetric N x N matrix (thus 
X^ — X and X^ — X). In the Hermitian case, the measure is given by: 

dX^YldXij . (4.47) 

i<j 

The partition function (4.46) and the measure (4.47) are invariant under the trans- 
formations: 

X^VXV'^ , (4.48) 

where V is an element of 5'0(iV, M). 

One way to see the aforementioned correspondence is by relating the loop equa- 
tions of the two models. The loop equation for the model (4.46) was extracted in 
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[18], and involves only the resolvent: 



uJz) = 4tr^-^ . (4.49) 
^' N z-X ^ ' 



This loop equation takes the form: 



ujs{zf-^uj'^{z)-2W\z)uJs{z) + 2~f,{z)) =0 , (4.50) 



where {■ ■ .)s denotes averages computed in the model (4.46) and fs{z) is the random 
polynomial: 

q W'(z) - W'(X) 

az):^^tr—^ . (4.51) 

^ ^ N z-X ^ ' 

Relation (4.50) corresponds to our equation (4.21) under the identifications: 

uj{z) < — > LJs{z) (4.52) 
f{z)^fs{z) . (4.53) 

Since the quantities k{z),m.Q{z), l{z) are determined by a solution of (4.21) via equa- 
tions (4.22), this gives a one to one correspondence between solutions of the two 
models' loop equations. 

A more direct relation between the two models can be extracted from their eigen- 
value representations. To see this, we must first discuss the eigenvalue representation 
of (4.46). 

In the context of the Dijkgraaf-Vafa correspondence, one must use the holomor- 
phic [4] version of (4.46). This is obtained by allowing X to be a complex symmetric 
matrix (thus removing the hermiticity constraint X^ — X <^=^ X — X) and consid- 
ering gauge transformations of the form (4.48), where now V is an invertible complex- 
valued matrix subject to the constraints V'^ = and detF = 1. This amounts to 
working with the complexified gauge group SO{N, C) (then the normalization prefac- 
tor IG^I is also modified as explained in [4]). In that case, (4.47) becomes the natural 
top holomorphic form on the complex space A/" = {X G Mat{N,C)\X'^ = X} and 
the integral in (4.46) must be performed along a real, boundary-less submanifold A 
of this space whose dimension equals the complex dimension of Af and whose closure 
is non-compact. The admissible choices of A are constrained by the requirement 
that A be stabilized by the action (4.48) of the complexified gauge group and that 
the integral (4.46) converge when calculated along A. This constrains the choice of 
A in terms of the leading coefficient of W [4] . 

To see the relation to (4.44) explicitly, it suffices to notice that integrating out 
the angular variables in (4.46) leads to: 

Zs= [d\i... [ dXj^H (A, - A,)e-T Sf=i m^i) ^ (4.54) 



i<3 
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where 7 is a suitable open contour in the complex plane (since 7 is constrained 
only by the leading term of W [4], it can be chosen to coincide with the contour 
used in (4.44)). The form (4.54) corresponds to choosing A := {X G J\f\ 3V G 
SO{N, C) such that VXV^ = diag(Ai . . . A^) with Ai . . . G 7}. Note that the 
real dimension of A equals the complex dimension of J\f, as required. Indeed, a 
generic complex symmetric matrix X can be diagonalized by a complex orthogonal 
transformation V G SO{N, C): 

VXV^ = diag(Ai . . . A^) (4.55) 

with complex Xj. Since this is true outside a complex codimension one locus in A/", 
it follows that imposing the constraints Xj G 7 simply halves the number of real 
parameters in X. 

The factor Y[i<j {■^i ~ ^j) is just the square root of the factor Y[i<j i^i ~ 
familiar from the case of the adjoint representation of GL{N , C). This can be seen 
most easily in the holomorphic matrix model set-up by repeating the argument given 
in Appendix A of [4] with the observation that the number of integration variables 
in (4.47) is reduced with respect to the case treated there due to the condition 
Xij — Xji. Obviously (4.54) coincides with the reduced matrix integral (4.44): 

Zred = Zg . (4.56) 
Using relation (4.42), we find that the partition function of (4.1) can be written as: 

Z = (27rz)*(*+^)/2^,ZQ , (4.57) 

where Zq is the quantity defined in equation (4.43). The factorization (4.57) shows 
that the average in our model of any functional -F($) which does not depend on S, 
A or Qf coincides with the average of F{X) in the SO{N) model (4.46): 

{Fm^{F{X)), . (4.58) 

This explains the identifications (4.52) and gives the precise relation between the two 
models. 

4.5 The resolvent loop equation from the eigenvalue representation 

To derive the loop equation of (4.44) (or, equivalently, of (4.46)), one can follow 
standard procedure by starting with the identity: 

N N 



1=1 1=1 L p<q 
Performing the partial derivative gives the relation: 



. (4.59) 



,A 1 N^W'(Xi) V- 1 1 ^ . . X 
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Introducing the traced resolvent (4.12) and using the identity: 

1 1 _ 1 1 1 V ^ 



allows us to write (4.60) in the form: 



(4.61) 



(ujIzY -uj Iz) ^ > — ) = 



N 



N ^ z - Xi 
1=1 



(4.62) 



Using the degree d — 1 polynomial (4.20): 



- g ^ W'(z) - ly'(Ai) 



leads to the loop equation in the form: 



{uizf - jt^'iz) - 2W'{z)cu{z) + 2f{z)) = 



(4.63) 



(4.64) 



which recovers (4.21). 



Consider the spectral density p(A) = i Xlili ^(-^ ~ '^i) expansions: 



j>0 



(4.65) 
(4.66) 
(4.67) 



In the large N limit, we have: 



uJo{z) :— lim {uj{z)) — g dX 

N^oo J 



Po(A) 



/o(^) := hm {f{z)) = g / dApo(A) 



W' {z)-W'{X) 
z-X 



AT— >oo 

In this limit, equation (4.62) becomes: 

uoizf - 2W'{z)u;o{z) + 2fo{z) = , 



(4.68) 



(4.69) 



which shows that the quantity ^0(2:) = i-Oo{z) — W'{z) is a branch of the hyperelliptic 
Riemann surface: 

u'^W\zf-2fo{z) . (4.70) 
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The loop equation (4.62) can also be written in the form: 

, / n9 9 // n f dx W'(x)u!(x) , „ /, „ X 

M.)^-4.'W-2y^-4^) = (4.71) 

where C is a counterclockwise contour encircling all of the eigenvalues of l» but 
not the point z. The form (4.71) can be used to give an iterative solution for the 

coefficients ujj{z) of the large expansion (4.65). In this paper, we are interested in 
the CP^ and RP^ diagram contributions u}o{z) and u}i{z). Expanding (4.71) to order 
0{g/N) and using the relation: 

{ojizr) = {u{z)f + 0{{g/Nr)=uo{zr + ^uJo{z)uj,{z)+0{^^^ , (4.72) 



we find: 



u;o(zr-2d>^^^^^^^^0 (4.73) 
^ ^ ^ 2m z - X ^ ^ 



c 

/V /V ^ f dx W'(x)uJi(x) ,, , /, 
2u;o{x)oui{x) -2d, — — ^-±-1^ = ^(x) . (4.74) 
J Q Zttz z — X 

The first relation is equivalent with the large N limit (4.69) of the loop equation 
and determines uJo{z) in terms of the polynomial fo{z). The second relation is an 
inhomogeneous integral equation for uJi{z), which constraints this quantity once /o(^) 
(and thus uJo{z)) has been fixed. 

5. Relation between the matrix model and field theory 

5.1 Compcirison of loop equations and Konishi constraints 

In this subsection we find a map from matrix model to field theory quantities which 
takes the loop equations (4.21) and (4.22) into the Konishi constraints (3.45), (3.46) 
and (3.47). We set Np — 2 and Np = 8 from now on. 

We start by considering a complex^ microcanonical ensemble for our holomorphic 
matrix model following the procedure of [4] . This is obtained by introducing complex 
chemical potentials Hk for a partition of the complex plane into domains with 
smooth boundary, followed by a Lcgendre transform of the resulting grand-canonical 

generating function J-'(t,fi) = — In Z{t,fi), which replaces the chemical po- 

tentials by complex variables Sk^^ (remember that tj/j are the coefficients of W). 
This produces the desired free energy (=microcanonical generating function) F{t, S), 
which satisfies the equations: 

OF , , 



case. 

10. 



Since our matrix model does not have a real counterpart, the set-up of [4] is essential in our 
These should not be confused with the eigenvalues of the symmetric matrix S! 
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Working with this microcanonical ensemble amounts to imposing the constraints: 



ifk) = S, (5.2) 
where fk is the fiUing fraction of the domain Dk, which we define through: 

A = £ |^^(.) . (5.3) 

Here is the boundary of D^. Note that the variables Sk must satisfy: 

J2s, = g , (5.4) 

k 

SO this a 'constrained' microcanonical ensemble, as discussed in more detail in [4]. 
Thus passing to the microcanonical ensemble allows us to eliminate g in terms of Si, 
if we choose to treat Si as independent variables. 
Expanding (5.2) to 0{g/N) inclusively gives: 

/ dz 

and: 

^^i(^) = 0. (5.6) 

Imposing these conditions fixes all coefficients of fo{z) and fi{z). In fact, con- 
dition (5.5) singles out one solution ujq of (4.73) while (5.6) specifies the associated 
solution of (4.74) by selecting the trivial solution of the associated linear homoge- 
neous equation. Then a simple counting argument along the lines of [21] shows that 
that {i^{z)) is completely determined as a function of tj and Sk- This also determines 
the expectation values of k{z),mQ{z) and l{z) via equations (4.22). 

Let us now consider the Konishi constraints (3.45) and (3.46) of the field theory, 
which determine all relevant quantities through equations (3.47). To specify uniquely 
a solution {R,T) of (3.45) and (3.46), we shall impose the constraints: 

Note that we must impose two conditions on the solutions of (3.45), (3.46), which 
amounts to fixing the coefficients of both the polynomials / and c. As in [7], the 
second equation in (5.7) can be viewed as a rigorous quantum definition of the rank of 
the k-ih. factor of the unbroken low energy gauge group. This interpretation requires 
that we choose such that they separate the critical points of W, which we shall 
assume from now on. 
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To map a solution of the matrix loop equations to a solution of the Konishi 
constraints, we shall identify: 

<Sk = Sk ■ (5.8) 

Using (5.4), we find that (5.8) fixes the value of the matrix coupling constant in 
terms of field theory data: 

g = J2^k (5.9) 

k 

Consider the large N expansions: 

{k{z)) = (^Ij m,{z), (5.11) 



i=o 

oo 



j=0 Vjv/ 

j=0 

Expanding (4.21) to leading order, we find: 

cooizf -2W'{z)u;o{z) + 2fo{z)^0 . (5.14) 
We also expand (4.22) to leading order and order g/N to obtain: 

rnqoiz) = -u'q{z) (5.15) 
lo{z) = u;o{z) . 

Comparing (5.14) with (3.45) shows that u!o{z) and fo{z) should be identified with 
R{z) and f{z) respectively. Moreover, equations (5.15) agree with the first and the 
last two equations in (3.47) provided that we identify K{z) with ko{z) as well as 
Mq{z) with 4mQo{z) and L{z) with Hq^z). 

To recover (3.46) and the second equation in (3.47), we consider the g/N terms 
of (4.21) and of the first equation in (4.22), which read: 

2uoiz)ui{z) - u'oiz) - 2W'{z)ui{z) + 2fi{z) = (5.16) 

and: ^ 

ki{z) ^ uq{z)ui{z) ^ -Uq{z) . (5.17) 
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Consider the operator: 

^:=E^*5| . (5.18) 

k 

Applying this to both sides of relation (5.14) and of the first equation in (5.15) gives: 
2uJo{z)Su;o{z) - 2W'{z)Suo{z) + 26fo{z) = (5.19) 

and: 

Sko{z) = uJoiz)5cuoiz) . (5.20) 
Combining these two equations with (5.16) and (5.17) leads to the relations: 

Wo(-2)[5c^o(-s)+4c^i(;^)]-l^'(;^)[^c^o(;^)+4c^i(;^)]-2c^o('S) + ['5/o('S)+4/i(;^)] = (5.21) 
and: 

Sko{z) + 4ki{z) = 2uJo{z) + ujo{z)[Sujo{z) + 4uJi{z)] . (5.22) 

These relations agree with (3.46) and (3.47) provided that wc identify T(z) with 
SiUQ{z) + AuJi{z) as well as M{z) with dkQ{z) + 4ki{z) and c{z) with Sfo{z) +4/i(z). 

In conclusion, matrix model and field theory quantities must be identified ac- 
cording to the table: 



Matrix Model 


Field Theory 


Sk 


Sk 


uoiz) 


R{z) 


6iUo{z) + Aui{z) 


T{z) 


6ko{z) + Aki{z) 


M{z) 


ko{z) 


K{z) 


AmQo{z) 


Mq{z) 


Ahiz) 


L{z) 


Hz) 


m 


Sfo{z)+iMz) 


c{z) 



Table 1: Identification between field theory and matrix model quantities. 

This correspondence recovers that used in [18] upon applying the field theory 
and matrix model relations (3.53) and (4.52), which connect our model with the 
SO{N) theory. 

5.2 The effective superpotential 

The identifications of the previous subsection allow us to determine the field theory 
effective superpotential up to a term independent of the coefficients tj of W. This 
contribution can be identified independently by using the Veneziano-Yankielowicz 
computation of Section 2.3. 
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For this, we note the relation (tr$-^) = jfr-, which imphes: 



H^)) = J^f (8.23) 



where ^pj^^^ Sj>o 'J+^w- Combing this with the expansion: 



^ = E(|)'''.- (5-24) 



gives: 



j>0 



On the other hand, one has the obvious field theory relation: 

(tr$^) (5.26) 



dtj j 
which gives: 

= WIT)''''" ■ ^'-''^ 

Using the identification T{z) — Su!o{z) + 4:U!i{z) (where 5 is the operator given 
in (5.18)), this implies: 



(5.28) 



Si — Si 



Here is a function which depends only Si but not on the coefficients of W^^. Since 
we always set Si — Si, we shall only use the notation Si from now on. Notice that we 
have derived (5.28) without having to postulate some analytic continuation which 
would avoid identifying integer or real quantities with complex numbers. This is 
because the use the formalism of [4], which automatically avoids such problems. Also 
notice the prefactor of 4 in front of the RP^ contribution Fi, which arises naturally 
in our derivation. The fact that diagrams of topology MP^ generally contribute with 
a factor 4 to the effective superpotential was previously discussed in [46] by using 
the perturbative superfield approach of [5]. 

Expression (5.28) determines the effective superpotential only up to the coupling- 
independent term il>{S). Together with the contribution to (5.28) from the non- 
perturbative part of F, this term should correspond to the Veneziano-Yankielowicz 
potential computed in Subsection 2.3, which cannot be determined through Konishi 



^^Remember that relation (5.9) eliminates g in terms of Si. 
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anomaly arguments. Applying the conjecture of [1, 47] to our model leads to the 
proposal: 

d 

^|JiS)^aJ2Si , (5.29) 

i=l 

where: 

a = (AT -4) In A (5.30) 

with A the field theory scale of Subsection 2.3. To check this expression, we now 
proceed to compute the non-perturbative contribution to W^ff in the matrix model 
and compare with the results of Subsection 2.3. This will allow us to complete the 
proof of the relation: 

dF 

^e// = 5]iV^^ + 4Fi + a5^5, . (5.31) 

i ' 1=1 

5.3 Computation of the Veneziano-Yankielowicz superpotential from the 
matrix model 

In this subsection we show how the Vcncziano-Yankiclowicz contribution to the ef- 
fective superpotential can be extracted from the matrix model. We shall follow the 
approach of [47, 7], by computing the non-perturbative contribution to the matrix 
integral and checking agreement between the non-perturbative part of (5.31) and the 
result (2.26) of Subsection 2.3. 

Let us consider the classical matrix vacuum: 

(l>) = diag(Ail^^,...,A,l^J , {S) = {A) = {Qf) = . (5.32) 

where Xj are the critical points of W. Following [47, 7] , we shall compute the Gaussian 
approximation to the matrix integral expanded around this vacuum. This is the 
semiclassical approximation in the background (5.32). 

Since we wish to compare with field theory, we must work in the microcanonical 
ensemble, which constraints g through relation (5.9). In the semiclassical approxi- 
mation about the background (5.32), one has Si = Ni by equations (5.2) and (5.3). 
Then (5.9) implies g = N, which means that the prefactor of the action in the ex- 
ponential of the matrix integrand must be set to one. In particular, the large N 
expansion can be reorganized in powers oil/N rather than g/ N. For simplicity, we 
can therefore start with: 

Z = e-^^ — ^ / d^dAdSdQe-'-^'^^^^^^^^^-^UQpQ^ (5.33) 
vol(G) J 

and impose the microcanonical ensemble conditions Si — after performing the 
semiclassical approximation. We have set Np — 2, since this is the only case when 
the matrix model is well-defined. 
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As in Subsections 4.3 and 4.4, the partition function (5.33) can be reduced to: 

Z = (27rz)^(^+i)/V^Z, , (5.34) 

where Zg is the partition function of the matrix model with SO{N, C) gauge group 
and a complex matrix X in the symmetric representation. 

It is clear from (5.34) that we can expand Z around this vacuum by expanding 
Zg around the background (X) = ($). We let x := X — (X) denote the fluctuations 
of X and decompose x into iVj x Nj blocks Xij. In the reduced model Zg, we have 
an SO{N, C) gauge symmetry which allows us to set the off-diagonal blocks of x to 
zero. Thus we can choose the gauge: 

Xij — for i ^ j . (5.35) 

To implement this in the BRST formalism, we introduce ghosts Cij, antighosts Cij 
and Lagrange multipliers Bij. These transform in the adjoint representation of 
SO{N,C), so Cfj = -Cji, Cjj = -Cji and Bfj = -Bji. The quadratic part of 
the gauge-fixing action is: 

d 

= XI i^^j^^^j + ^i^i - ^j)CijCij) , (5.36) 

where we used the BRST transformations: 

sX^i[C,X], sC^iC^, sC^B, sB^O . (5.37) 

When expanding to second order in x, the diagonal blocks Xi := Xa acquire masses 
rrii = W"{Xi). Hence the non-perturbative piece of the partition function (5.33) is 
given by: 

Here G' = nti SO{Ni) is the unbroken gauge group. Since the action is quadratic, 
we can choose the integration contour: 

BijeiR xaeR, (5.39) 

which gives: 

^np = n[(27r)^'^^(A.-A,)^»^^]n 

i<j i=l 



(27ri) 



NjiNj + l) 



'n 



Ni 



yol{SO{Ni)) \W"{X^ 



TT 



(.5.40) 
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Using [20]: 



iV^ N N 2N 
log[vo\{SO{N))] = -— log ^-^ + _ log — + 0(iV°) , (5.41) 



we extract the free energy Fnp — — log Znp'. 



np 



E 

i=l 



... (-47r3) N? 
— H loK log 



— H 

4: 4: 



Ni N,. , 2N, 



27re3/2 '4 tt^ 



-J2NiNj\og[27r{\-Xj)] 



(5.42) 



As explained above, we have Si = Ni since we are in the semiclassical approximation 
of the matrix model. Replacing by Si in (5.42), we can identify the quadratic 
terms in Si with F^^ and the hnear terms with Ff^: 



S? 



i=l 
d 



s,. 



i=l 



' {-2ti^)S, 
_e3/2pF"(Ai) 

SiW"{\i) 



- J]5,cS, log[27r(A,-A,)] , (5.43) 

(5.44) 



According to (5.31), the non-perturbative contribution to the effective superpotential 
should be given by: 



1=1 * i=l 



(5.45) 



Ignoring constant and linear terms in Si (which can be absorbed by a finite renor- 
malization of the field theory scale A of Subsection 2.3), we find: 



1=1 



■A'''-'w"{Xi)'^'^'lliXi-Xjy 



-2Ni 



s 



Ni-2 



+ 0{S) , (5.46) 



which recovers the leading contribution (2.26) derived by threshold matching. 



6. Conclusions 

We studied a class of non-anomalous, chiral M = 1 U (N) gauge theories with an- 
tisymmetric, conjugate symmetric, adjoint and fundamental matter in the context 
of the Dijkgraaf-Vafa correspondence. By using the method of generalized Konishi 
anomalies, we extracted a set of chiral ring constraints which allowed us to identify 
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the 'dual' holomorphic matrix model and give an explicit expression for the gaugino 
superpotential in terms of matrix model data. This gives a proof of the Dijkgraaf-Vafa 
correspondence for a class of theories with quite nontrivial chiral matter content. 

As a by-product of this analysis, we found that the effective superpotential of 
our models coincides with that produced upon confinement in non-chiral SO{N) field 
theories with a single chiral superfield transforming in the symmetric representation. 
This provides an independent proof of a relation suggested by holomorphy arguments. 

Our results encourage us to think that similar methods could be applied success- 
fully in order to extract non-perturbative information about more realistic models 
employed in supersymmetric phenomenology. It would be interesting to see how far 
this program can be implemented. 

Surprisingly, we found that the number of fundamental flavors Np in field theory 
must be taken to differ from the number of flavors Np in the dual matrix model. This 
is unlike the non-chiral case considered in [11, 13, 14, 16, 15, 39, 8], for which the 
number of fundamental flavors agrees between field theory and the matrix model. 
More precisely, consistency of the matrix model requires Np — 2, while anomaly 
cancellation in our field theory requires Np — 8. Despite this disagreement, one 
can match the Konishi constraints in the chiral ring with the loop equations of the 
two-fiavor matrix model. 

Another interesting question concerns the geometric engineering of our models, 
which provides their embedding in JIB string theory and leads to an interpretation 
of the gaugino superpotential as the flux-oricntifold superpotential of certain non- 
compact Calabi-Yau backgrounds. As in the case of the U (N) theory with symmetric 
or antisymmetric matter [21, 22, 24] (whose geometric engineering was discussed in 
[25] ) , this can be achieved by applying the methods of [37] . A deeper understanding of 
the geometric realization could also shed light on the mismatch between the number 
of flavors in the matrix model and fleld theory. A detailed investigation of these 
issues is under way [36]. 
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A. Gauge invariance of the matrix model measure 

With an arbitrary number of flavors Np, we have the following transformations under 
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the action (4.7): 

dA {detuf~^dA 

dS ^ {detU)-^^+^^dS (A.l) 
dQ (detU)^^dQ , 

where U is an arbitrary element of GL[N, C). Therefore, the matrix model measure 
dfjL — d^dAdSdQ transforms as: 

d/i {detuf^-'^dii , (A.2) 

and is invariant if and only if Np = 2. 

The first relation in (A.l) is familiar from the (adjoint) holomorphic one- matrix 
model [4], while the last relation is obvious. Thus wc only have to prove the second 
and third equation. For this, it suffices to check them for diagonalizable U, since our 
matter representation is continuous and the set of diagonalizable elements is dense 
in GL{N, C). Thus we can take U = VTV-\ where V is an element of GL{N, C) 
and T — diag(ii . . .^jy), with Yij^j 0- The explicit form of the transformations 
(4.7) gives: 

dA [det Ra{U)] dA 

dS [det Rcs{U)] d^ , (A.3) 

where Ra and Res are the antisymmetric and conjugate symmetric representations 
of GL{N,C): 

Ra{U){A) = UAU^ , A^ =-A 

R,s{U){S) ^U-^SU-^ ^U-^S{U-^Y , S'^^S . (A.4) 

Using Ra{VTV-^) = Ra{V)Ra{T)Ra{V)-^ and R,s{VTV-^) = Rcs{V)Rcs{T)R,s{V)-\ 
we find: 

dA [det Ra{T)] dA 

dS [det Rcs{T)] dS . (A.5) 
Since T is diagonal, one easily obtains: 

deti?„(r)= H ktj = {Ylt,f-' 

det RUT) = n ^^"^•' = (11^^)"^'*^'^ ' 

SO that: 

det Ra{U) = {detUf-^ and det Rcs{U) = (detC/)"^"^ . (A.7) 
This leads immediately to the second and third equation in (A.l). 
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B. Classical vacua of the SO{N) model 



For the SO{N) model, the only matter is the complex superfield X, subject to the 
condition X^ — X and transforming as follows under the gauge group: 

X VXV^ = VXV-^ . (B.l) 

Here V is valued in SO{N,M.). The F-flatness relations for Wtree = triy(X) read: 

W'(X) = (B.2) 

while the D-flatness condition has the form: 

1[X,X]=0 , (B.3) 

where the left hand side is the moment map^^ for the action (B.l) of SO{N,'R) 
(computed with respect to the natural symplectic form U!{X, Y) — — tr {XY) on the 
representation space). 

Writing X = + iX*™, where X^'^ and X'"* are two real- valued symmetric 
matrices, the D-flatness constraint becomes: 

[j^re^j^im] , (B.4) 

which shows that X*^'^ and X*"* can be diagonalized simultaneously by performing a 
gauge transformation (B.l): 

X- = diag(AriiVi...AIflAr,) 

X^- = diag(AfliVi...A^^ljv,) . (B.5) 

Here the real numbers A^^ and A^'" are such that the pairs (A^'^, A^"*) are mutually 
distinct. 

Writing A^ := A^.^ + iA^™, equations (B.6) show that X can be diagonalized via 
a gauge transformation (B.l): 

X = diag(AiliVi ■ ■ ■ XdInd) , (B.6) 

where A^ arc distinct complex numbers. 

The F-flatncss condition (B.2) now shows that we can take D = d = degW'{z) 
and Xi . . . \d must coincide with the distinct critical points of W (again we use the 
convention that the block Afcl^Vj. does not appear in (B.6) if = 0). In such a 
vacuum, the gauge group is broken down to [nLi OiNk)]/{-l, 1}. 

As usual, we used dualization with respect to the KiUing form rj) = — tr (^77) of o{N) in order 
to view the moment map as being valued in the Lie algebra o{N) rather than in its dual. 
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C. Derivation of the generalized Konishi constraints 



Let us outline the derivation of the generahzed Konishi relations. We will make 
heavy use of the chiral ring relations (3.5) and (3.7)-(3.10). 

The first two equations (3.22-3.23) are straightforward modifications of the re- 
lations found in [7] for the U{N) model with a single adjoint field $ and no extra 
matter. The only differences are due to the coupling of $ to the symmetric and 
antisymmetric fields and appear on the left hand sides of these relations. 

The derivation of the other relations is somewhat more involved. Let us consider 
first the symmetric field S and the general transformation: 

SS - X^SY (C.l) 

where X and Y are arbitrary Grassmann even or odd N x N matrices. We start by 
computing the anomaly term induced by this transformation. For this, recall that 
Wq acts on 5S according to the conjugate symmetric representation: 

Wa.SS ^ -W'^SS - 5SWa , (C.2) 

where juxtaposition on the right hand side stands for matrix multiplication. 
Applying this formula twice, we find the anomaly term: 

-327rM := Tr (>V".>V..^) = 

Qgij -r yVas ^ ^YV r q^^^ yVa,s + YV r yVa,m g^^j 

= l-ti {X)tr (YW^) + (XYW^) + tr (W°X)tr (YWa) + 

Zi 

(-)[^ltr (XW^FWa) + -tr (W2X)tr (F) + ^tr (XWV) . (C.3) 

where \X\ denotes the Grassmann parity of X, Tr (. . . ) stands for the trace in the 
conjugate symmetric representation and tr (. . . ) is the trace in the fundamental rep- 
resentation. 

Applying this relation for X = Y = ^3^, we obtain: 



327r2 



1 ^ W2 Wo 
tr — tr — + 2tr — r + tr — tr 



$ z-^ {z-^Y z-^ z-^ 



= T(z)R(z) - 2R'(^) - lw"w„ . (G.4) 

The left hand side of the generalized Konishi relation (3.12) involves the variation 
of the super potential contracted with ^S: 

■ dW 1/1 $^\l/$ 1 \ 
^^''^ = -otr —s]+-tr -A —S] + (G.5) 

+gj^^^^^g/ = tr [s^-A] ^Q'^^S^-Qf . 
^z-^^ z-^ V z-^ J ^z-^^ z-^ 
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To arrive at this equation, we used the identity (4.30). Combining (C.5) and (C.4) 
leads to the Konishi relation (3.26). 

We next consider the transformation: 



5 A = XAY'^ . 



Since: 



Wc,.6A = WJA + 6AyVi , 
this leads to the anomaly term: 

85 A\ 



(C.6) 
(C.7) 



-327v'A := Tr >V".>V«. 



dA 

dAii dAj-i 



dA, ' 



^tr (X)tr (YW^) - (XYW^) + tr (W"X)tr (YWa) - 
{-f^\x (XWFWa) + ^tr (W2X)tr (F) - ^tr (XW^F) . 



(C.8) 



Choosing X — Y = one finds 

1 



327r2 



1 ^ 
tr tr 2tr ■ 



tr tr 



z - ^ z-^ 



= T{z)R{z) + 2R'{z) - -w"Wa . 



(C.9) 



On the other hand, the choice X = and Y — gives: 



-32tt^A = tr 



z-^ J 

where we used the spinor identities: 



(C.IO) 



(C.ll) 



The left hand side of the anomaly equation for X = Y = has the form: 



dW 



OA 



tr S 



z - ^ 



-A 



(C.12) 



where again we used identity (4.30). Together with (C.9). this leads to the Konishi 
relation (3.24). 
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If we chose X — and Y — , the left hand side becomes: 



We next notice that: 

I S^A^] . -U (J^A^Q^SA ^ -t. (S^A^ I ^ 

(C.14) 

where in the first hne we used invariance of the trace under transposition and in the 
second hne we used the symmetry of S and antisymmetry of A together with the 
spinor identity V'^V'/^ = —ipfjip^. A similar derivation shows that: 

I S^Am.\ ^ t. is^A- 



2;-$ z-^'^ j I z - ^ 2;-$^ 

where in the last line we used the chiral ring relation AW^ = —WpA. Combining 
this with the anomaly term (CIO) leads to the Konishi relation (3.25). 

Let us now turn to transformations of the fundamental fields. Starting with: 

sQf = J2^f^v^Q9 , (C.16) 



z-^ 

9=1 



we find: 

Np 



2 ^ Xf.Qpj^Q, ^ B.{z)Xff (C.17) 

5=1 

where we do not sum over the flavor index /. This is the most general form of the 
Konishi anomaly without any restriction on Xfg. Since Xfg is arbitrary, we conclude: 

'^Qpj^Q, = ^(z)Sfg . (CAS) 

Taking the trace of this equation leads to the Konishi relation (3.28). 
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